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We derive an effective low-energy Hamiltonian for potassium loaded zeolite A, a unique ferromagnet 
from non-magnetic elements. We perform ab initio density functional calculations and construct maximally 
localized Wannier functions for low-energy states made from potassium s electrons. The resulting Wannier 
orbitals, spreading widely in the alminosilicate cage, are found to be the superatomic s and p orbitals in 
the confining potential formed by the host cage. We then make a tight-binding model for these superatomic 
orbitals and introduce interaction parameters such as the Hubbard U. After mean-field calculations for 
the effective model, we find that ab initio spin density functional results are well reproduced by choosing 
appropriate sets of the interaction parameters. The interaction parameters turn out to be as large as the 
band width, ~ 0.5 eV, indicating the importance of electron correlation, and that the present system is an 
interesting analog of correlated multi-orbital transition metal oxides. 
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1. Introduction 

Electron correlation in materials, which often causes 
various non-trivial phenomena, has been one of the cen- 
tral issues in condensed matter physics. While there are 
many kinds of fascinating strongly correlated systems 
such as organic molecular p-, transition metal d-, or 
heavy fermion /-electron systems, it is of great interest to 
consider the possibility of correlated s-electron systems. 
Since valence s-electrons are usually itinerant, we expect 
that correlations in s-electron systems cannot be so sig- 
nificant. However, recently, alkali-metal loaded zeolite is 
attracting broad interests as a unique exception; when 
clusters of alkali atoms are confined in the cage of almi- 
nosilicate, their kinetic energy is drastically suppressed 
and the system can exhibit a variety of correlation effects. 

Although the unitccll of alkali-metal loaded zeolites 
is huge and complicated [typically it contains 0(100) 
atoms], the low-energy electronic structure of zeolites is 
surprisingly simple. The host alminosilicate cage makes 
a gap of several eV, in which the guest alkali s-electrons 
form a few narrow bands around the Fermi level. This 
situation is described by the so-called superatom pic- 
ture, 1-3 ) where each zeolite cage is regarded huge 
atom. Here, the confining potential formed by the host 
framework and the guest alkali-metal s electrons cor- 
respond to the atomic potential and valence electrons 
of the superatom, respectively. Recent ab initio calcu- 
lations have elucidated that the band dispersion around 
the Fermi level is well represented by a very simple tight- 
binding model, 3 ~ 8 ) suggesting that correlation effects in 
low-energy physics of zeolites can be systematically de- 
scribed by the many-body superatom model. 

As for ab initio derivation of many-body superatom 
models for zeolites, there have been several attempts 
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for alkali-metal-loaded sodalite. 4-8 * 1 Sodalite is one of the 
most prototypical zeolites made from the smallest build- 
ing block of the alminosilicate cage (the so-called /3-cage) . 
Each cage has only one valence s electron, so that it can 
be viewed as a crystal of hydrogen-like superatoms. In 
ref. 4, a single-orbital extended Hubbard model for the 
superatomic s electrons has been constructed, and the 
correlation strength was found to be much larger than 
the band width. 

On the other hand, we also have many intriguing su- 
peratomic multi-orbital systems in the family of zeolites. 
Potassium loaded zeolite A (Ki 6 Sii2Ali 2 048, hereafter 
we call it K-LTA, where LTA is the code for the crystal 
structure of zeolite A) is a typical example, for which 
a spin-polarized ground state is realized. 1,2 - 1 While the 
mechanism of spin-polarization has not been fully un- 
derstood, there are several experimental 9-12 ) and theo- 
retical 3, 13 ) indications that orbital degrees of freedom 
of the superatomic p orbital is important. In this pa- 
per, we derive an effective multi-orbital Hamiltonian of 
K-LTA from first principles. We estimate onsite poten- 
tials, transfer integrals, and interaction parameters in the 
Hamiltonian. We show that the energy scale of the elec- 
tron correlation, estimated as ~ 0.5eV, is as large as the 
band width, so that the system is strongly correlated and 
can be regarded as a unique "miniature" of transition 
metal oxides. 

The structure of the paper is as follows. In § 2, we first 
briefly describe the atomic configuration of K-LTA and 
give the level diagram of the low-energy states, which is 
useful to discuss the band structure of K-LTA in terms 
of the superatom picture. In § 3, we describe the detail 
of the method to evaluate the parameters in the effec- 
tive Hamiltonian. The basic idea is following: We first 
perform density functional calculation and make max- 
imally localized Wannier functions 14 - 1 for several bands 
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around the Fermi level. Using these Wannier functions, 
we then construct a tight-binding model and introduce 
interaction parameters such as the Hubbard U. Choosing 
appropriate sets of the interaction parameters, we per- 
form a mean-field calculation for the effective model and 
fit to the result of ab initio spin density functional calcu- 
lation. In § 4, we present our results which indicate that 
the system is indeed classified strongly correlated 
system. Section 5 is devoted to a summary and outlook. 

2. Level Diagram of Superatomic p-orbitals in 
Zeolite A 

Let us first look at the atomic configuration of K-LTA. 
The system has two kinds of alminosihcate cage, the so- 
called a and (3 cage (see Fig. 1 and related descriptions 
in ref. 13). While the diameter of the (3 cage is ~ 7 A, 
the a cage is much larger (~ 11 A), so that the guest 
potassium cluster is accommodated only in the a cage. 
In the unitcell, there are two inequivalent a and (3 cages, 
which arc denoted as ot\, ai, Pi, and fa- Concerning the 
atomic positions of potassium atoms, neutron scattering 
measurements 15 ) revealed that there are four kinds of 
sites: Kj, K n , K m , and K IV [Fig. 1 (a)]. Here, the K I; K n , 
and Km atoms sit around the center of the six-, eight-, 
and four-membered rings in the alminosihcate network, 
respectively, and the Kiv atom sits at the center of the 
a cage. The sites of Kt and Kn are fully occupied, while 
the Km and Kpv sites are partially occupied. Following 
the argument in ref. 13, we assume that the system has 
three-fold symmetry. 

The low energy bands around the Fermi level are 
mainly made from the superatomic p orbitals in the a 
cages. When we discuss the superatomic p levels, there 
are two key parameters. One is the difference between 
the potentials of the two a cages, defined as 



A7=min{e PBai ,ep^ ai ,e Pj , ai }— mia{e Pa , a2 ,e 



}(1) 



with e p ai specifying the ionization potential of the su- 
peratomic p^ orbital in the cage. Since the a\ cage ac- 
commodates more K + ions than the ai cage, we naively 
expect that the potential of a\ is deeper than that of ai. 
However, it should be noted here that the Ki atoms [de- 
noted by solid (blue) circles in Fig. 1 (b)] are distributed 
sparsely in the a.\ cage than in the ai cage, and the ol\ 
cage does not have Kiy at the cage center [solid light 
gray (pink) circles in Fig. 1 (c)]. These facts make, on 
the contrary to the above expectation, the ai-cage po- 
tential shallower than that of ai [the schematic diagram 
in Fig. 1 (b)]. 

The other important parameter is the crystal field 
splitting which is controlled by the atomic configuration 
of Km [solid dark gray (green) circles in Fig. 1 (c)]. For 
the Km configuration spreading in the xy plane [(i) and 
(iii) in the figure] , the degenerated p x and p y levels be- 
come lower (lower schematic diagrams). In contrast, for 
the configuration which extends in the z direction [(ii) 
and (iv)] , the crystal field lowers the p z level. Note that 
the ai cage contains 6 Km atoms and the ai cage does 
3 Km atoms. Also note that the Kiv atom is accom- 
modated only in the ai cage. The crystal-field splitting 



AV ai in the a, cage is defined as 



(2) 
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Fig. 1. (Color online) (a) Atomic configuration of K-LTA, (b) 
Difference in the cage-potentials (A/) formed by Kj clusters [de- 
noted by solid (blue) circles], and (c) Crystal-field splitting (AV ai ) 
formed by clusters consisting of Km and Kiv [denoted by solid 
dark gray (green) and light gray (pink) circles, respectively]. The 
z axis is taken to be perpendicular to one of the six-membered 
rings in the a cage, and we assume that the system has three- 
fold symmetry along the z axis. Note that the sign and size of the 
crystal-field splitting depends on the atomic configuration of the 
cluster. Four types of the Km and Krv configurations considered 
in the present study are depicted. 



In the previous study, 13 -* we performed spin density 
functional calculations for the four geometries in Fig. 2, 
which we call geometries I- IV. The figure also displays 
the schematic level diagrams with the superatomic pic- 
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ture. It was found that the geometry I, whose total en- 
ergy is the lowest among the four geometries, has the 
largest magnetic moment (~ 2/^) in the «2 cage. On 
the other hand, the magnetic moment is smallest for the 
geometry II. This geometry dependence was simply dis- 
cussed in terms of AI and {AV^} in ref. 13. While these 
are key parameters for the formation of local spin in 
the a cage, transfers between the neighboring a cages 
is also important to discuss the magnetism in the bulk. 
We should also consider the degrees of freedom of the 
/3 cage neglected in the previous analysis, mediating the 
hopping between the superatomic states in the a cage. 
On top of that, careful analysis on the interplay between 
electron itinerancy and correlations are imperative. Thus 
in the present study, we evaluate the interaction param- 
eters such as Hubbard U for the superatomic orbitals. 

3. Effective Model Construction 

Let us move on to the detail of the method to derive 
an effective model from first principles. For evaluation 
of interaction parameters such as the Hubbard U , in the 
previous study for alkali-metal loaded sodalite, 4 ) we used 
the constrained random phase approximation. 16, 17 ) How- 
ever, since the system size of K-LTA is much larger than 
that of sodalite and the numerical cost is so expensive, 
here we take a different approach. 

First, we introduce the following model Hamiltonian, 



H = Ho + Hint , 

h q =J2 e 



(3) 



^int — ^ ^ Uif-fiTlifji^Tliui^ ~\- ^ ^ U , ,i rtirn 'f^i7ncr'^'im / a 



i . n i 

E 

i : m<m f ,a 



>)n 



truer 'W (7 J 



(5) 



where c\ ma {ci miy ) creates (annihilates) an electron with 
spin <7 in orbital m at site i, and rii m(J = c| mo .Ci mcr . The 
term Ho describes the kinetic-energy part of the effective 
model, where £i m and U m j rn i are the ionization potential 
and transfer, respectively. In this model, we take into ac- 
count of not only the superatom p states in the a cage 
but also the superatom ,s states in the j3 cage, which 
lie just above the Fermi level as explained in the next 
section. The term H m t is the interaction part, i.e., the 
intraorbital (Ui m ) and interorbital {U' imm ,) Coulomb in- 
teractions and the Ising component of Hund's coupling 
{Jimm')- Since the standard spin density functional cal- 
culation is not SU(2) symmetric, we ignore the spin-flip 
component of the Hund's coupling and the pair-hopping 
interaction. For simplicity, we also ignore orbital depen- 
dence of interaction parameters in the a cage; common 
U a is used for any m and common U' a and J a are used 
for any combination of m and m' . For the /? cage, there 
is no orbital degree of freedom. 

Now, let us consider to derive the interaction parame- 
ters {U a , Up, U' a , and J a }. These parameters are deter- 
mined by the fitting procedure in such that the Hartree- 
like mean field solution of the model reproduces both of 



ab initio non-magnetic and spin-polarized density func- 
tional results based on the generalized gradient approxi- 
mation (GGA). 

In order to relate the mean field in the model calcula- 
tion and the density functional calculation, we introduce 
the following two terms. One is a term to scale the self 
interaction (SI) and the other is a term which guarantees 
that the non-magnetic model calculation gives the same 
result as the non-magnetic density functional calculation. 
For the first issue, as is well known, GGA includes the SI 
effect. 18 -' On the other hand, the Hubbard model defined 
in eq. (3) does not have SI; the intra-orbital interaction 
between electrons with the same spin are excluded in the 
model Hamiltonian. Therefore, to mimic the GGA situ- 
ation in the model calculation, we need to introduce the 
SI term as 



n S i[p] = J2 S ^ 



(0) 



where p lmana ' = {c\ na ,c ima ) is the density matrix and 
the parameter Si scales the size of SI. It should be noted 
that %si[/°] is a functional of the diagonal part of the 
density matrix. With this treatment, we can estimate 
the size of SI in GGA by %si[pgga], where pgga is the 
density matrix obtained from the ab initio non-magnetic 
GGA calculations. For the second issue, we introduce 
71[pgga\, where "H[p\ denotes the Hartrec approximation 
for Hi n t , being a functional of the density matrix p. The 
"H[pgga] term is evaluated with the density obtained by 
non-magnetic GGA. 

Thus the one-body part Ho in cq. (4) is defined as 

Ho = Hgga - Hsi[pgga] - H[pgga], 

where Hqga is the tight-binding representation of the 
non-magnetic GGA Kohn-Sham Hamiltonian with e^ A 
and t^fin' in diagonal and off-diagonal parts, respec- 
tively. We use the maximally localized Wannier functions 
for the basis functions for those matrix elements. More 
precisely, the matrix elements of Ho are calculated as 



,GGA 



Ui(Si + 1) 



.GGA 



GGA\ 
r i rn rn ) 



(7) 



and 



GGA 



Hmjm' — timjm' • (8) 

Here, pf^n 1S the matrix element of the non-magnetic 
GGA density matrix and iVf GA = J2 m P?,nm ■ Notc tha t 
transfers in Ho is the same as those in Hgga, since Hsi 
and H are local in real space. 

Now, the mean- field Hamiltonian %mf[p] for H reads 

as 



Hmf[p}=Ho + H[p}+Hsi[p}- 



(9) 



Note that the result of non-magnetic GGA is exactly 
reproduced as Hmf[pgga] in the present framework. 
We then consider to solve "Hmf by spin-polarized self- 
consistent Hartrec scheme. We can estimate the intcrac- 



Q 



J. Phys. Soc. Jpn. 



FULL PAPERS 



Geometry I Geometry II 

ot] a 2 aj a 2 a x a 2 a] a 2 




Geometry III Geometry IV 

ot-i a 2 0C1 a 2 a 2 aj a 2 




Fig. 2. (Color online) Geometries considered in the present study and their expected schematic level diagrams of the superatomic-p 
orbitals. 
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tion and self-interaction parameters, U a , Up, U' a , J a , S a , 
and Sp, by fitting the model band dispersion to the ab 
initio spin density functional results. 



4. Results 

Using the Tokyo Ab initio Program Package, 19 ' we first 
performed density functional calculations with the GGA 
exchange-correlation functional. 20 ^ We used a plane- wave 
basis set with the ultrasoft pseudopotcntials. 21 ) The en- 
ergy cutoff was set to 36 Ry for wavefunction and 144 
Ry for spin and charge densities and 4x4x4 k-mesh was 
adopted. Maximally localized Wannier function with the 
ultrasoft pseudopotential was constructed following the 
recipe of ref. 22. 

4-1 Maximally localized Wannier functions 

The GGA band structures for the geometries ITV are 
shown in Fig. 3 as (red) dotted lines. We see that the 
low-energy band structures depend sensitively on the ge- 
ometry. Especially, the bandwidth is rather different from 
each other. From the six (seven) bands around the Fermi 
level, we constructed the maximally localized Wannier 
functions for the geometry I and II (III and IV), which 
are employed as bases of the effective Hamiltonian in 
eq. (3). We estimated tight-binding parameters for the 
geometries I-IV by calculating matrix elements of Hgga 
in the Wannier basis. In Fig. 3, we plot interpolated band 
dispersions [(green) solid lines], from which we see that 
original ab initio bands are well reproduced by the Wan- 
nier interpolation. The band around —0.4^—0.3 eV is 
formed by the superatomic s orbital in the a cage, which 
is not included in the effective model. 

As shown in detail later (Fig. 6), the s orbital in the a 
cage is fully occupied and its exchange splitting is negli- 
gibly small. Thus, for simplicity, we do not consider ex- 
plicitly this degree of freedom in the low-energy model. It 
should be noted that the s state in the a cage has small 
transfer integrals between the surrounding sites and tiny 
onsite hybridizations with the p states, because it is well 
localized in the cage and the confining potential is almost 
spherical. Moreover, the s state should have larger U and 
S than the p states, so that there should be a significant 
difference in the ionization potential, £i m [eq. (7)] , for the 
s and p states. Therefore, the s state is expected to be 
irrelevant for the magnetic properties. 

In Fig. 4, we display the isosurfacc contours of the 
plane- wave part of the Wannier functions for the geome- 
try I. In this geometry, there are six bases (s/?i, sj32,p z ot\, 
Pz&2, PyOt2, and p x cx.2), and the p x a\ and p y a\ states are 
not included in the effective model, because the latter 
two exist far above the Fermi level by, at least, ^0.5 eV. 
Similarly, for the geometries II-IV, the high-energy p^ai 
states were not included in the effective models. Table I 
lists our calculated spatial spreads of the Wannier func- 
tions, f2 m j, for all the geometries. As is indicated in the 
table, the size of the Wannier functions is as large as the 
diameter of the host cage (~11 A for the a cage and ^7 
A for the P cage) and flp < fl ma , so that the onsite in- 



teractions in the /3 cages should be larger than those in 
the a cages. 




Fig. 4. (Color online) Isosurfacc contours of maximally localized 
Wannier functions for supcratomic-p (-s) orbitals in the a (/3) cage 
of the geometry I. The amplitudes of the contour surface are 0.015 
[light gray (yellow)] and —0.015 [dark gray (blue)] in the atomic 
unit. Dark gray (green) and light gray (pink) spheres represent 
potassium Km and Kiv, respectively. 



Table I. Spatial spread of maximally localized Wannier func- 
tions for superatomic states, Q, and sclf-consistcntly determined 
ionization potentials, e defined by eq. (7). The difference in cage 
potentials AI [eq. (1)] and the crystal-field splitting in the ai cage 
AV a2 [eq. (2)] are also listed. The unit is given as A for Q and eV 
for e, AI, and AV a2 . 
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4-2 Effective low- energy Hamiltonian 

Let us move on to the details of the effective model 
derived with the present scheme. First, the transfer in- 
tegrals can be read from the table of the matrix el- 
ements of Hgga in the Wannier basis [eq. (8)]. The 
nearest-neighbor hoppings between the a and /3 cages are 
\tp* /vas p\ = 0.01 - 0.03eV and \t^ asf} \ = 0.02 ~ 0.04eV. 
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Fig. 3. (Color online) GGA band dispersions [(red) dotted lines] and interpolated band dispersions for the maximally localized Wannier 
functions [(green) solid lines] for geometry I-IV. Energy is referenced to the Fermi level. 
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On the other hand, the transfers between neighboring 
a cages are \t P:s0tPz0l \ = 0.01 ~ 0.02eV, \t PzaPa:/y a\ = 
0.02 - 0.03eV, and \t Px/y0iPxly0 \ = 0.03 - 0.04eV. The 
distant transfers are negligibly small. We remark that, 
among the transfers between a cages, the p z -p z hopping 
is smaller than those of other p states. This is because 
the p z orbital is not pointing at the neighboring a cages. 
Thus, the p z states tend to be localized and to have large 
magnetic moments. 

In order to estimate {U a , Up, U' a , J a , S a , Sp}, we 
performed spin-polarized self-consistent calculations for 
Hmf in eq. (9), and looked for the parameter set which 
reproduces the band dispersion of ab initio spin density 
functional calculations. First, we searched for the pa- 
rameter set for which the difference between the GGA 
band dispersion of the geometry I and that of the model 
calculation along the symmetry line W— >X— >■ V — >L— > 
K in the Brillouin zone is minimized. We assumed that 
U a < Up < 0.6eV, S a < Sp < 1, U'JU a < 1, and 
JajU a < 0.2 with the realistic parameter range. 4 ) To 
understand and explain this system by a minimum set 
of parameters, we then introduced a scaling factor x for 
{U a , Up, U' a , J a } and minimize the band-dispersion dif- 
ference in GGA and model for the geometry II, III, and 
IV. We found that the GGA band structures for these 
three geometries are successfully reproduced for x = 0.6, 
Sp = 0.7, and S a = 0.1. These parameters show that 
Geom. II-IV have smaller interaction parameters than 
Geom. I. This is consistent with the fact that the band 
widths in Geom. II-IV are broader than that in Geom. I 
as shown in Fig. 3. 

Before presenting the resulting interaction parameters, 
we summarize in Table I our calculated self-consistent 
non-magnetic ionization potentials, e in cq. (7). Using 
the e data, we estimated the effective cage-potential dif- 
ferences AI in eq. (1) and the crystal field splitting in 
the «2 cage AV a2 in eq. (2), which are given in the bot- 
tom two rows in Table I. 23 ) It should be noted here that 
these one-body quantities have appreciable geometrical 
dependence, suggesting that slight changes in the potas- 
sium configurations can cause substantial changes in the 
low-energy electronic structure of K-LTA. According to 
Table I, we draw in Fig. 5 the level diagram of the su- 
peratomic states in the a and (3 cages. We can see that 
the Q2-cage potential is significantly deep compared to 
the other cage potentials and thus will form the trapping 
potentials for the low-energy electrons. 

We next show in Fig. 6 our calculated spin-polarized 
model band dispersion of the majority spin [(red) solid 
lines] and minority spin [(green) dashed lines], compar- 
ing with those of the ab initio results [majority (+) and 
minority (x)]. Table II compares the resulting local mag- 
netic moments in the model and ab initio calculations. In 
spite of the limited number of the free parameters, the ab 
initio band dispersion and moment size are well repro- 
duced. 13 ' The magnetic moments are basically formed in 
the a 2 cage, being consistent with the level diagram in 
Fig. 5. The transfers cause the antiferromagnetic corre- 
lation between the a\ and ai cages. 

The interaction parameters estimated by the present 
fitting are summarized in Table III. We see that the re- 



sulting U is 0.24-0.5 cV and is always about ten times 
larger than the transfer integrals as 0.01-0.04 eV, thus in- 
dicating that the system resides in moderately /strongly 
correlation regime. We note that interaction parameters 
depend sensitively on the atomic configuration; they tend 
to be large (small) for the geometry I (II-IV). 24 ' It is 
also interesting to note that the self-interaction param- 
eter Si tends to be larger when the interaction param- 
eters are large. This is a reasonable trend, 18 - 1 because 
self-interaction is generally large for localized states, and 
localized states usually have strong onsite electron cor- 
relation. In the present study, the Si value is estimated 
by parameter fitting to the spin polarized GGA calcula- 
tions. Therefore, the present Si parameter qualitatively 
measures the localization of the spin-polarized wavefunc- 
tions, but not that of the paramagnetic wavefunctions. 



Table II. Local magnetic moment (M, in the unit of ^b) in the 
ab initio and the model calculation for K-LTA. 
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Table III. Estimated interaction parameters (U, U' , J) and self- 
interaction parameter (S) in K-LTA 
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Finally, let us comment on the controversial situation 
in the experiment. Experimentally, two kinds of possible 
magnetic structures have been proposed for ferromag- 
netism of K-LTA. One is the ferrimagnctic model, 9 ' and 
the other is the spin-canted antiferromagnetic model. 11 ' 
In the former, the system consists of two sublattices with 
large (~ 2.8 /^b) and small (almost 0.0 /j,b) magnetic mo- 
ment, while, in the latter, every a cage has a moment of 
1 and the spins are canted in the low temperature 
by the spin-orbit interaction. These two magnetic struc- 
tures seem to agree with the results for the geometries I 
and IV, respectively, 13 ' although the size of the moments 
is relatively smaller than the experimental values. It is 
interesting to note that the geometries I and IV are ener- 
getically stable (slightly lower in the geometry I by ~0.1 
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Fig. 5. Level diagrams for the supcratomic p levels in the a cages and the s levels in the f3 cages for the geometries I-IV. Energy zero 
is the Fermi level (dotted line). The difference between the cage potentials of a\ and 0:2, A J in eq. (1), and the crystal- field splitting in 
eti cage, AVa 2 in eq. (2), are also depicted in the geometry I. 
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Fig. 6. (Color online) Band dispersions of majority [(red) +] and minority [(green) x] spins obtained by the spin-polarized GGA, and 
those of majority [(red) solid line] and minority spins [(green) dashed line] calculated by the model calculation for the geometries I-IV. 
Energy is reference to the Fermi level. 
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e V). 13 ) The present results seem to support the former 
proposal, but we should emphasize that K-LTA resides in 
moderately or strongly correlation regime so that care- 
ful analyses with proper consideration of the correlation 
effect are required to clarify the mechanism of the spin- 
polarization in real K-LTA. Indeed, strongly correlated 
multi-orbital systems generally have rich phase diagrams. 
For example, a subtle competition between the interac- 
tion parameters and the tight-binding parameters can 
drastically change the magnetic properties of the system. 
Elaborated many-body analyses for K-LTA are indeed a 
fascinating future problem. 

5. Summary and Outlook 

We have derived an effective low-energy Hamilto- 
nian for K-LTA. First, we constructed a tight-binding 
model from GGA using the maximally localized Wan- 
nier functions. We then introduced the interaction and 
self-interaction parameters, and solved the model self- 
consistcntly. With appropriate sets of these parameters, 
we have succeeded in reproducing the result of spin den- 
sity functional calculation. We have found that the sys- 
tem reside in the strong coupling regime in that the 
bandwidth and electron correlation have the same en- 
ergy scale. We have also shown that the details of the 
potassium-cluster configurations crucially affects both 
the one-body and two-body parts of the effective Hamil- 
tonian of K-LTA. 

While the present method of evaluating the interac- 
tion parameters has a great advantage in that it can 
be applied to huge systems such as K-LTA, it does not 
work for systems for which spin density functional calcu- 
lation gives a non-magnetic ground state. Since the size 
of magnetic moment tends to be underestimated in GGA, 
the present method is expected to give the lower bound 
of the interaction parameters. It is indeed an interest- 
ing future challenge to apply other ab initio methods to 
evaluate interaction parameters such as constrained den- 
sity functional theory (ref. 25) or constrained random- 
phase approximation (refs. 16 and 17) to K-LTA. While 
formidable numerical cost will be required, comprehen- 
sive model construction by different methods is impor- 
tant to establish the low-energy correlation physics in the 
zeolite system. 
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